Limitations of light delay and storage times in EIT experiments with 

condensates 
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We investigate the limitations arising from atomic collisions on the storage and delay 
times of probe pulses in EIT experiments. We find that the atomic collisions can be 
described by an effective decay rate that limits storage and delay times. We calculate the 
momentum and temperature dependence of the decay rate and find that it is necessary 
to excite atoms at a particular momentum depending on temperature and spacing of the 
energy levels involved in order to minimize the decoherence effects of atomic collisions. 
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Electromagnetically induced transparency 
(EIT) [JIJ allows a medium that is opaque to be 
made transparent around its resonant frequency. 
Recent experiments using EIT of cold trapped 
atomic gases have shown that one can delay light 
pulses by reducing the group velocity to speeds as 
low as 1 m/s ||k|. Moreover it is possible to store 
them for up to 1 ms before regenerating them H 
(for a recent review see Q). There are two main 
factors that affect the photon delay or storage time 
in EIT experiments pg . Firstly, the delay time is 
limited by the finite frequency width of this trans- 
parency window. We shall assume that the pulse 
duration is sufficiently long that its bandwidth fits 
within the transparency window. Secondly, the de- 
cay of the dark state arising from atomic collisions, 
as well as from possible couplings to states outside 
the system, limits the delay and storage times in 
these experiments. In this paper we discuss how 
the decay of the dark state is related to the delay 
or storage times. We study the dependence of the 
decay width introduced by collisional atomic in- 
teractions on temperature and on the momentum 
transferred from the photons to the atoms. 

We consider a homogeneous gas of atoms whose 
internal structure may be described as a three-level 
system as shown in Figure Homogeneity should 
be a good approximation for the central region 
of the large alkali-metal condensates achieved in 
present day experiments. The ground state \B, 0), 
which describes a coherently populated, weakly in- 
teracting Bose-Einstein condensate with a mean 
number of N atoms, is resonantly coupled to the 
excited state \A) by a weak probe pulse. This 
probe pulse is described by a single-mode quan- 
tized electromagnetic field with frequency u) v . Its 
dipolar coupling strength to the atoms is g. For 
convenience, we take the probe pulse to be ini- 
tially in a number state \p; n) with n photons. The 
metastable translational state |C, k) of atoms with 
cm. momentum k in internal state \C) is reso- 




nantly coupled to \A) with a classical laser field 
with frequency to c characterized by the Rabi fre- 
quency ft. We describe the decay of the excited 
state | A) due to spontaneous emission with a rate 
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FIG. 1. We consider 3-level atoms illuminated by 
a single quantized (probe) mode with coupling g and a 
classical laser field characterized by the Rabi frequency 
f2. We assume the condensate atoms to be in a coher- 
ent state \B, 0). The lasers transfer the atoms to the 
internal level C, into the translational state |C, k) with 
defined momentum k. 7 a and jc(k) characterize the 
decay rates from |^4) and \C,k), respectively. 

The population of the state \C,k) decays pre- 
dominantly due to 2-body elastic collisions with 
the atoms in internal state \B) which conserve the 
internal atomic quantum number. This decay may 
be described by a rate 7c(&) depending on the mo- 
mentum of the decaying mode if the modes with 
momenta different from k as well as the condensate 
of atoms in internal level \B) are close to thermal 
equilibrium. jc(k) may be calculated perturba- 
tively using a linear response approach to the mode 
operator time evolution . 

The dynamics of the system we consider is de- 
scribed by the following system of equations of mo- 
tion for the atomic operators A, Ck and the probe 
photon operator a: 

A = —iflCk — igV No a, 

6 k = -inA-{ lc (k)/2)Ck-Fc,k, 
'a = -ig*fW A-{' 1A /2)a-F A . (1) 



Here we have chosen an interaction picture where 
the fast evolution due to the energy splitting of the 
states has been factored out of the operators, and 



the resonance conditions 6ab — £cb — ? 
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€ab — hio P have been taken into account. Here, 
ej? = h 2 k 2 /(2m) is the cm. kinetic energy of an 
atom in mode k and internal state \C). J^c.fe and 
Fa are the Langevin noise operators correspond- 
ing to the decay of the operators Ck and a re- 
spectively which provide for the mode operators 
to obey Bose commutation relations and the cor- 
relation functions of these operators to decay to 
their equilibrium values ^,0). 

The recently achieved light propagation slowing 
and storage of light pulses in cold atomic gases 
rely on the possibility to create a dark state. A 
dark state \D(n,8)) which is consistent with our 
assumption of a probe beam being initially in a 
number state and a condensate in a coherent state 
\B, 0) has the form of a "beam-splitter state": 



\D(n,e)) = (cos 
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where |C, k; m) is the number state where m atoms 
are in mode k in level C. If the mixing angle is de- 
fined through 



tan# 



(3) 



and provided one neglects jc(k), the state (||) does 
not change in time except for an overall phase. By 
adiabatically adjusting Q and thereby varying 9 
from to tt/2, one can transfer the initial pho- 
ton population into collective atomic excitations 
in \C,k). The reverse can be achieved by adiabati- 
cally changing 9 back to 0. In |||| this is described 
in terms of dark-state polaritons. In this ideal- 
ized process the sum of photons in the probe beam 
state and atoms in \C,k), i.e. n p (t) + nc,k(t) — 
(a 1 (i)a(t)) + {Cl(t)Ck(t)) , remains constant and 
the model describes the stopping and regenerating 
of a light pulse in an atomic medium. 

This simple picture no longer fully holds once 
decay is taken into account. We have studied 
the effects of Jc(k) by numerically calculating 
Ti p (t) + nc,k{t)- To this end we use the Heisen- 
berg equations of motion (]]]) with respect to an 
initial dark state (||), assuming the correspond- 
ing energies to have small imaginary components 
as given in ([!]). Note that we have in princi- 
ple to take into consideration the corresponding 
Langevin noise operators Fc,ki Fa when calculat- 
ing the time evolution of the number operators. 



One may show, however, that a quantum regres- 
sion theorem holds for the correlation functions of 
two normal ordered operators 0. Therefore ([!]) 
may be used, neglecting the Langevin terms, for 
the calculation of their time evolution. 

We see in Figure |^ that n p + nc,k decays faster 
the closer 9 is to tt/2, or equivalently, SI is to zero. 
This faster decay arises because more photons are 
converted into atomic excitations. Apart from this, 
the decay is dominated by jc(k) and as long as this 
is reasonably small (i.e. Jc(k) < 9"\/Nq> which is 
well satisfied in present experiments, e.g. ■yc(k) < 
(27r)10kHz, 5 V^o - (2tt)15MHz in §j), the pop- 
ulations rip + nc,k, n p and nc,k all decay at the 
same rate. 
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FIG. 2. Number of photons in the weak probe 
beam plus the number of atomic excitations in \C,k), 
n p + nc,k, as a function of 6 for different times us- 
ing experimental parameters Q|: Q. = (27r)5.61 MHz, 
7A = (27r)10 MHz. We also assume 7c(fc) < gs/Na. 
(N ~ 10 6 atoms, n ~ 3 • 10 4 photons.) 



In the light storage experiments, the system re- 
mains at 9 = tt/2 (i.e. fi = 0) for much longer 
than it spends in transition. For example, in Q 
the delay time plus the time it took to turn on and 
off the coupling laser was on the order of microsec- 
onds whereas the maximum storage time was on 
the order of milliseconds. Since the rate of decay is 
highest when the photons are stored (9 = tt/2) and 
the storage time far exceeds the delay and turn- 
on/off times, we approximate that all of the decay 
of n p + nc t k occurs while the photons are stored as 
atomic excitations. The maximum storage time t s 
is then given by t s < l/jc(k). 

For slow light experiments, the characteristic 
decay time, Td, at which n p + nc,k has decayed 
to a fraction 1/e of its initial value n, has to 
be calculated numerically. Assuming that effects 
other than collisions are irrelevant |lC(| , we esti- 
mate Td for coherent slow light propagation in a 
Bose-condensed gas at T = from the time evolu- 



tion of the dark state (|2j) , for the experimental pa- 
rameters from |[: Q = (2tt)5.61 MHz and gV^o - 
(2tt)10MHz. We find that r d ps l.2/j c (k). As can 
be seen in Figure ^ Td is in general greater than 
l/ 7 c(fc). 

We have seen that the decay caused by colli- 
sions and characterized by 7c (fc) limits the max- 
imum time for the storage of light pulses in an 
atomic condensate. We have calculated 7c (k) 
for the decay of a single excited mode of atoms 
in a different internal state than the condensate 
atoms, using stationary perturbation theory in the 
interactions beyond the quadratic Hartree-Fock- 
Bogoliubov hamiltonian. The results of this calcu- 
lation can also be derived using a time-dependent 
linear response approach [Q. In the following we 
study the dependence of "fc(k), and hence of and 
r s , on momentum k and temperature. We consider 
a homogeneous dilute condensate of Nq particles 
dominated by two-body collisions. Their strength 
is determined by the s-wave scattering length a, 
which we assume to be the same for collisions ei- 
ther of atoms within level B or between atoms in 
B and C . An extension of our methods to the case 
of trapped condensates is straightforward. We as- 
sume the condensate is confined in a box of volume 
V and apply periodic boundary conditions. Colli- 
sions between atoms in \B) change the free par- 
ticle dispersion relation to the Bogoliubov quasi- 
particle dispersion relation 

~4 = ^y k ^2 + y 2 k (4) 

where /1 — (hko) 2 /2m and yt = k/k^. The in- 
verse of the healing length, fc = {^Tm^a) 1 / 2 , sep- 
arates the phonon-type from the free-particle type 
region of the quasi-particle spectrum. To the lead- 
ing order, only the spin-conserving collisions be- 
tween atoms in \B) and |C) affect the excitations 
in |C, k), causing a non-zero 7c (k). Collisions in- 
volving atoms in \A) can be ignored because of its 
negligible population. Finally, collisions between 
atoms in internal state |C) can be neglected as 
there is no macroscopically occupied cm. mode. 
We neglect the energy shift of the levels due to 
interactions between excitations in |C) and \B), 
which arises at the same order as 'yc(k) @- 

The collisional decay rate can be separated into 
the Beliaev-like and Landau- like decay rates (here- 
after referred to as Beliaev and Landau decay 
rates), i.e. jc(k) = 7c( fc ) + 7c( fc )- Tne Beliaev 
decay rate, corresponding to a particle in |C, k) col- 
liding with a condensate atom in \B, 0), producing 
one quasiparticle in \B,i) and one particle excita- 
tion in \C,j), both of lower energy, zf — zf + zf 

{zf = y^2 + y 2 , zf = y 2 k ), is given by 



Jo \ + {zff ) 

x(l + nf + nf). (5) 

with populations nf — (exp[zf n/ksT] — 
and nf = (exp[(zf + z C b - zf)^/k B T] - l) -1 . 
(zcb = £cb/ 70- The Landau decay rate is caused 
by the reverse processes where a particle in |C, k) 
collides with a quasiparticle in \B,j) to produce 
a particle excitation in |C, j) with higher energy 
zf = zf + zf. It is given by 

7&(*)=^r^fi— r ^=) 

x(nf-nf), (6) 

where nf = (exp[{zf + z CB + zf)^/k B T] - I)- 1 . 
Whereas Beliaev processes are restricted to modes 
whose energies are less than that of decaying parti- 
cles, Landau processes can involve any mode with 
an appreciable population. This causes the Lan- 
dau decay to dominate at higher temperatures. At 
T = only the Beliaev rate survives giving 

m t (i ln (^ + v /T +^A 

7c(«0t=o = auj k y k 1 3 ■ 

V Vk J 

(7) 

At high momentum where > fcp, 7c(fc)T=o ~ 
"fKT(k) 7 where ^KTfk) = uqctv is the collision rate 
from classical kinetic theory, assuming the parti- 
cles have speed v = Uk /m and cross sectional area 
a = Ana 2 . At low momentum k <C fc , we obtain 
7c(^)t=o ~ hk° /96mirnQ. This low-momentum 
limit of the decay rate has the same k 5 depen- 
dence as the T = Beliaev decay rate for quasi- 
particle excitations in a single state condensate 
[11 13], from which it only differs by a constant 
numerical factor. 

In Figure ^| we plot 'jc(k) as a function of k/ko 
for different temperatures using experimental pa- 
rameters of Q . As the temperature increases from 
T = 0, Landau processes start to dominate at 
low k while Beliaev processes remain dominant at 
high fc, creating a minimum decay rate at a finite 
k. Eventually the minimum disappears because 
at that point Landau processes dominate at all k. 
The minimum of 7c (fc) at a finite k develops be- 
cause of the separation ecB of the energy levels 
\B,0) and |C, 0). Note that for the condensate in 
Q with a length of 230 /zm the minimum y k = k/ko 
is of the order of 5 • 10 -3 , such that the description 



of the loss from \C, k) by a decay rate is expected 
to be justified. 
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FIG. 3. The decay rate 7c (k) due to atomic 
collisions, as a function of normalized momentum 
k/ko for different temperatures, using the experi- 
mental parameters M: 23 Na atoms, a = 2.8 nm, 



T c = 435 nK, n 



10" cm" 3 , fi/h = (2tt) 1 .21 kHz 



and ecis/h = (27r)f.8GHz. 

To emphasize the importance of ecs, we plot 
7p(fc), Jc(k), and 7c (k) as functions of k for differ- 
ent values of ecB in Figure |[ We see that, as tcB 
is increased, the Landau rate significantly increases 
at low k because a large ecB causes (nf — nf ) 
in eqn. (||) to approach while (pM — nf) ap- 
proaches zero as k — > for — > 0. Further- 
more, the Beliaev rate decreases slightly as ecs 
increases because nf w in eqn. (§) at large 
We can see that the minimum in 7c develops for 
sufficiently large ecB at certain temperatures be- 
cause the Landau rate dominates at low k while the 
Beliaev rate dominates at high k. As the tempera- 
ture increases, the momentum transfer required to 
minimize 7c (k) also increases. For example, using 
T/T c = 0.5 and the experimental parameters used 
in Figure |j, we find t s — 0.5 ms at y k — 0.1 and 
t s = 4.1 ms at y k = 1. 

We have shown in this paper that the storage 
time for light pulses in atomic condensates using 
dark-state based EIT schemes is limited by loss 
caused by collisions between particles and quasi- 
particles in the atomic ensemble. We found that 
one can minimize this loss in a homogeneous BEC 
by adjusting the momentum k transferred to the 
atoms. The optimum k depends on the characteris- 
tics of the three-level system and the temperature. 
In this respect the collisional decay in a system 
where the atoms in the excited translational mode 
are in a different internal state than the conden- 
sate atoms shows an important property distinct 
from the case with only one internal state. 
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FIG. 4. The decay rate yc(k), as a function of 
k/ko f° r different values of zqb = ecs/M, f° r the 
same experimental parameters as used in Figure |^ and 
T/T c = 0.1. We chose z CB = 0.0001 (dash-dotted 
line), zqb = 0.01 (dashed line), and zcb > 10, where 
~ (solid line). The two insets show "/c(k) and 
7c (k) for the same parameters. 
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